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r9 ■ Abstract 
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Renormalization factors for three- and four-quark operators, which appear 
in the low energy effective Lagrangian of the proton decay and the weak in- 
teractions, are perturbatively calculated in domain-wall QCD. We find that 
the operators are multiplicatively renormalizable up to one-loop level without 
mixing with any other operators that have different chiral structures. As an 

q^ , application, we evaluate a renormalization factor for Bk at the parameters 

where previous simulations have been performed, and find one-loop correc- 

"tt \ tions to Bk are 1-5% in these cases. 
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I. INTRODUCTION 

Calculation of hadron matrix elements of phenomenological interest represents an in- 
evitable application of lattice QCD. In the past decade much efforts have been devoted for 
the calculation of three- and four-quark hadron matrix elements relevant to the proton de- 
cay amplitude and the weak interaction ones using the Wilson and the Kogut-Susskind(KS) 
quark actions. However, the satisfactorily precise measurement of the matrix elements has 
not been achieved so far because of the inherent defects in these quark actions: the explicit 
chiral symmetry breaking in the Wilson quark action causes the non-trivial operator mixing 
between different chiralities and for the KS quark it is hard to treat the heavy-light cases 
due to the flavor symmetry breaking. 

The domain- wall quark formulation in lattice QCD, which is based on the introduction 
of many heavy regulator fields, was proposed by Shamir |l]|l anticipating superior features 
over other quark formulations: no need of the fine tuning to realize the chiral limit and no 
restriction for the number of flavors. Recent simulation results seem to support the former 
feature non-perturbatively 0-§|. It is also perturbatively shown that the massless mode at 
the tree level still remains stable against the quantum correction 0. These advantageous 
features fascinate us to the application of the domain-wall quark for calculation of the three- 
and four-quark hadron matrix elements. 

In order to convert the matrix elements obtained by lattice simulations to those defined 
in some continuum renormalization scheme (e.g., MS), we must know the renormalization 
factors connecting the lattice operators to the continuum counterparts defined in some 
renormalization scheme. In this article we make a perturbative calculation of the renor- 
malization factors for the three- and four-quark operators consisting of physical quark fields 
in the domain- wall QCD(DWQCD). This work is an extension of the previous paper |7j], in 
which we developed a perturbative renormalization procedure for DWQCD demonstrating 
the calculation of the renormalization factors for quark wave function, mass and bilinear 
operators. We focus on whether or not the renormalization of the three- and four-quark op- 
erators in DWQCD is free from the notorious operator mixing problem. In the Wilson case 
it is well known that the mixing problem is not adequately manipulated by the perturbation 
theory, leading to an "incorrect" value for the Bk matrix element. 

This paper is organized as follows. In Sec. [TT] we briefly introduce the DWQCD action 
and the Feynman rules relevant for the present calculation to make this paper self-contained. 
In Sec. pl| our calculational procedure of the renormalization factors for the four-quark 
operators is described in detail. We also evaluate the renormalization factors for three- 



quark operators in Sec. IV. In Sees. Ill and IV numerical results for one-loop coefficients 



of the renormalization factors are given with and without the mean field improvement. In 
Sec. [V], using our results, we analyze a renormalization factor for Bk- Our conclusions are 
summarized in Sec. [VI]. 

The physical quantities are expressed in lattice units and the lattice spacing a is sup- 
pressed unless necessary. We take SU(iV) gauge group with the gauge coupling g and the 

N 2 -l 
second Casimir Cp = , while iV = 3 is specified in the numerical calculations. 

2N ' P 



II. ACTION AND FEYNMAN RULES 

We take the Shamir's domain- wall fermion action [l|], 

- J2 (i>(n) s (-r + jf,)Uf,(n)ip(n + n) s + ip(n) s (-r - i^U^n - n)ip{n - /i) s ) 



N 3 

Sdw = Yl Yl 

n s=l 



2 , 



+- (V»(n) a (l + j 5 )ip(n) s+ i + tp(n) s (l - 75 )V(n) s _i) + (M - 1 + 4r)^(n)W(n)« 
+ m^ ('0(n)j Vs P + '0(w)i + ip(n) 1 P^ip(n) N: 

n 

where n is a four dimensional space-time coordinate and s is an extra fifth dimensional or 
"flavor" index, the Dirac "mass" M is a parameter of the theory which we set < M < 2 
to realize the massless fermion at tree level, misa physical quark mass, and the Wilson 
parameter is set to r = — 1. It is important to notice that we have boundaries for the 
flavor space; 1 < s < N s . In our one-loop calculation we will take N s — > oo limit to avoid 
complications arising from the finite N s . Pr/l is a projection matrix Pr/l = (1 ±7s)/2. For 
the gauge part we employ a standard four dimensional Wilson plaquette action and assume 
no gauge interaction along the fifth dimension. 

In the DWQCD the zero mode of domain-wall fermion is extracted by the "physical" 
quark field defined by the boundary fermions 

q{n) = P R i>(n) 1 + P L ijj(n) Ns , 

q(n)=^(n) Ns P R + ^(n) 1 P L . (2) 

We will consider the QCD operators constructed from this quark fields, since this field has 
been actually used in the previous simulations. Moreover our renormalization procedure is 
based on the Green functions consisting of only the "physical" quark fields, in which we 
have found that the renormalization becomes simple ]7|]. 

Weak coupling perturbation theory is developed by expanding the action in terms of 
gauge coupling. The gluon propagator can be written as 

G A J(k) = M ^_ 3 ^__ (3) 

in the Feynman gauge with the infrared cut-off A 2 , where sin 2 (fc/2) = J] M sin 2 (/c M /2). Quark- 
gluon vertices are also identical to those in the N s flavor Wilson fermion. We need only one 
gluon vertex for our present calculation: 

V*(k,p) 8t = -tgT A {^cosi-k^/2 + P J2) - irsm(-kj2 + pj2)}5 sU (4) 

where k and p represent incoming momentum into the vertex (see Fig. 1 of Ref. 0). T A 
(A— 1, . . . , iV 2 — 1) is a generator of color SU(iV). 

The fermion propagator originally takes N s x N s matrix form in s-flavor space. In the 
present one-loop calculation, however, we do not need the whole matrix elements because 
we consider Green functions consisting of the physical quark fields. The relevant fermion 
propagators are restricted to following three types: 



(q(-p)q(p)) 



-i^fn sinpn + (1 — We a ) m 



Sqip), 



-(1 -e a W)+m 2 (l -We- a 
(q(-pW(p, s)) = ^ (i llx sinp, - m (l - We'*)) (e^^P* + e^ 8 " 1 ^) 



+ F 



m li'-f^ sin p^ — m f 1 — We 



e~ a (V q(s_1) 



P R + e 



-a(N s -s) 



Pi 



(il>(-p, s)q(p)) = j (e-^ N ^P L + e-^'-VPn) (^ sinp, - m (l - We" Q )) 



(5) 



(6) 



(7) 



with 



W = l-M-r^(l-cosp M ), 



cosh(o; 



l + W^ 2 + E M sm 2 p M 



2|W| 
F = l-e Q W-ro 2 (l - HV 



(8) 

(9) 
(10) 

where the argument p in the factors a and W is suppressed. 

In the perturbative calculation of Green functions the external quark momenta and 
masses are assumed to be much smaller than the lattice cut-off, so that we can expand the 
external quark propagators in terms of them. We have the following expressions as leading 
term of the expansion: 



l-wl 



(qq)(p) = - 

if + (1 — WQJm 

(# s )(p) = (qq)(p) (w^Pl + wZ'-Pr) 
m){p) = (w^Pr + wZ'-Pl) (qq)(p), 



(11) 
(12) 
(13) 



where wq = 1 — M. 



III. RENORMALIZATION FACTORS FOR FOUR-QUARK OPERATORS 

We consider the following four-quark operators: 

1 r 



O, 



(qil^i.qzl^qi) ± (gi7„ 94) (£57^2) 



A„,L„ \(- rpA„,R, 



1 = -CAqaMiqslM + (qiT A lM(QsT A lM, 



1 
2N' 



°2 = 7^7(?i7^ 2 )(?37» + (gi^7^2)(g 3 T A 7^ 4 ), 



(14) 
(15) 
(16) 



where 7^ ,R = ^^Pl^r- Summation over repeated indices such as /i and A is assumed. We 
note that q,i (i = 1,2,3,4) are boundary quark fields in DWQCD. For the convenience of 
calculation we rewrite the above operators as 



(17) 
(18) 
(19) 

where a, b, c, d are color indices, and 0, represent the tensor structures in the color space: 

ab;cd 



o ± = 
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ab;cd 


XftTZdm-fidX 
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ti)(q 
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SnhS, 



ab:cd 



ab^cdi 



SadScb- 



(20) 
(21) 



To derive these formula, we have used the Fierz transformation for 0± and the formula 



E^ 



1 


r l _ ~ 1 




- — 101 + 101 


2 


iV 



(22) 



for 1)2 . 

We calculate the following Green function: 

\ij\kl 



(Or)Xs = (WMtoYfiWm) 



(23) 



where T = ±, 1, 2. Spinor indices are labeled by a, (3, 7, 5 and color ones by i, j, k, I. Truncat- 
ing the external quark propagators from (Or), where we multiply (Or) by ifa + (1 — wl)m, 
we obtain the vertex functions, which is written in the following form up to the one-loop 
level 



(l-» T(Ar)^ = (l-»„Y(A!? , +A< 1) )^ 



y./3 ;7<5 



(24) 



where the superscript (i) refers to the i-th loop level and the trivial factor (1— Wq) a is factored 
out for the convenience. We suppress the external momenta pi since the renormalization 
factor does not depend on them. 

The tree level vertex functions A r are given by 



r = ±, 
r = i, 
r = 2, 



1 
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7j®7j" 


a/3;7<5 . 


1 

2 


'tf®l*~ 


a/3;7<5 



101 ±101 



ij;kl 



1 r 



%®rf 



a/3;7<5 



-iVl0l + l0l 

ij:kl 

101 ' 



ij;kl 



(25) 
(26) 
(27) 



where acts on the Dirac spinor space representing [jx 7y]a/3;7<5 = (7x)a/3(7y)75- 

The one-loop vertex corrections are illustrated by six diagrams in Fig. 1, the sum of 
which yields the one-loop level vertex function 



Af> 



<■* d 4 k 



/_ " (/?+ ,..., +/ ,). 



(28) 



In order to obtain the expressions for the integrands Jp, . . . , Jp we should note that the 
internal quark propagators appearing in the diagrams are multiplied by the damping factor 



which comes from eqs. (|12|) and (|13D. The following formula are useful 
(# s > (w'q^Pl + w^- s P R ) = (w' - x P R + w^~ s P L ) faq) 






(# s > (w'^Pr + w^~ s P L ) = (w s Q - l P L + vtf'- a P R ) m) 






(29) 
(30) 



where F = e a — W and F = e a — Wq. Here we set m = Pi = for the internal propagator. 
The contribution from Fig. la takes the form 

G = \ja B {v,(k)G(k) + V,(k)G(k)}r x [G(k)V u (k) + G{k)V v {k)} <g> T Y G^(k), (31) 



where Tx = 7^, IV = 7^ or 7^, and the interaction vertices are 
y ii = -igifn cos(k^/2), Vp = -rg sin(/c M /2). 



(32) 



The color factors are represented by J^ B , which are listed in Table |I[ In a similar way the 
contributions from Fig. lb and Fig. lc are given by 

^r = \jb B {v,(k)G(k) + V,(k)G(k)}T x ® [V v (-k)G(-k) + V v {-k)G{-k)}v Y G^(k), (33) 






\jc B {v^k)G{k) + V^k)G(k)}T x ® T y {G(k)V u (k) + G(k)K(k)} G^(k). 



After a little algebra the expressions of /£' ' c are reduced to 
1 



Ir = ^9 2 Ja A K [T + A VA ] [T x 



Ty] 

,2/ 



*r = -7,9 JbK [TT X ® T Y + cos\kj2) sin 2 k a ( lilla T x ) <g> (t^IY) 



Iv = ^9 J< K [TT X ®T Y + cos\kj2) sin 2 k a ( lflla T x ) 



(r y 7 Q 7 M ) 



where 



K 



F 2 F 2 (Asm 2 (k/2) + X 2 ) 
T = r 2 sm 2 (k/2)F 2 + r sin 2 kF, 
Ava = Yl cos 2 (k fl /2) sin 2 fc M . 



(34) 



(35) 
(36) 
(37) 

(38) 

(39) 
(40) 



In order to rewrite the second term of 7 r ' we apply the Fierz transformation 
(7^7*7^) ® ilvlctlu) = -lu ®lu=l L ® it, 

(lulali) ® ilvlalu) = -1 L V © 7^(1 - 2<J QI/ )(1 - 2<V) 



-1 L V ® 7, L (1 - 2<W(1 - 2^) + 2 7, L ® 1 L J. 



afij 



(7m7«7^) ® (7M7a7?) = 2P R © P ^m« = lu ® 7*<W 
(7/x7a7* ) ® (7^7«7m) = 2P/? ®P L = lv®ll 



& 



(41) 
(42) 
(43) 
(44) 
(45) 



where [jx 7y]a/3;7<5 = ( r Yx)as('lY)'y/3, and summation over v is taken. The Fierz transfor- 
mation is again used for the second equality. We omit the tensor term in eq. flHD since it 
vanishes in the integral. 

Choosing Y X y = 7^ in eqs. (|35"D , (|36"D and ([37]) we first consider the case of 0±. After 
simplifying the expressions of the color factors we obtain 



4 = -g 2 K(T + A VA ) 



lu®lu 



( ± 



l ± 



--g 2 K(T + A S p) 



lu®lu 



-g 2 K{T + A 



VA) 



ii®ii 



(C F ±i)l®l T 2^101 
.^Li®l + (I ±(7^)101 



(46) 

(47) 
(48) 



The total contribution becomes 

1 



11 + 11 + n 



g 2 K 



lu®ll 



af3\~j8 



101 + 101 



ij;kl 



x{TC F + A VA {C F + ±-±) 



. ,1 1 ,-. 

Asp[ 2~2N )l 



(49) 



and 



ja +I b 



-g 2 K 



lv lu 



a/3;-fS 



x {TC F + A VA {C F 



101-101 
1 



ij\kl 



2\' +Asp{ 2 + 2N )} > 



(50) 



where 



A SP = cos 2 (/c/2) • sin 2 k. 



We should note that the other three contributions 7" ' ' c from Fig. la', 16' and lc' are equal 
to Jp' ,c respectively, therefore the factors 1/2 in eqs. ( 
contributions of all. 

Comparing the one-loop results to the tree level ones we obtain 



(51) 

ig. la', lb' and lc' are e 
and (E(ft disappear in the total 



Ah 
A. 



1 + 9 
1 + 9 



,N 



N 
2 N + 1 



N 



{< T > (N + 1)+ < A VA >(N + 2)- < A SP >} 
{< T > (N - 1)+ < A VA >(N- 2)+ < A SP >} 



A (o) 

A 



+ > 

(0) 



with 



<X>- 



d A k 

-7r(2^ 



K{k)X{k) 



(52) 
(53) 

(54) 



for X = T, A VA , A S p. We remark that C F < T + A VA > and C F < T + A S p > correspond to 
the one-loop vertex corrections to the (axial) vector current and the (pseudo) scalar density 
which are expressed as (T VA — 1) and (T S p — 1) in Ref. JFj. The expressions of eqs.(|52"D 
and ([5*3"D show an important property of 0± in the DWQCD formalism: the one-loop vertex 



corrections are multiplicative. This is contrary to the Wilson case, in which the mixing 
operators with different chiralities appears at the one- loop level [§]. 

We next turn to the case of 0\^- For (D\ the vertex corrections of eqs.(|35]), (|36|) and (p7|) 
with T x = 7^ and T Y = 7^ are written as 

1 



11 = -g 2 K(T + A VA ) 
l\ = - l - g *K(T + A VA ) 
I{ = l -g 2 K{T + A SP ) 



lv®lv 



lu ® lu 



lu®lu 



'-Nc F + l -)mi-^\Qi 

A lx ~ , N 1 . ~ 
Il®l + (-^ + Cp.)l0l 



(55) 
(56) 
(57) 



The total contribution including those from Fig. la', lb' and Id is given by 



2{I a 1 +I b 1 + I< i )=2-g 2 K 



x 






a/9;75 



J Up -|- A.V a — — Acp 

v 2 2N 



-M<8>1 + 1©1 
1 



ij;kl 



(58) 



Using the tree level result in eq. (p6|) the vertex function up to the one-loop level is expressed 
as 



Ai 



1 



l + g 2 -{<T> (N 2 - 1)+ < A v A > N 2 - < A SP >} 



Ai 



(0) 



(59) 



This result shows that the operator 0\ is multiplicatively renormalizable in DWQCD, which 
is in contrast with the Wilson case ||. 

In a similar way we write the vertex corrections for Oi- 



I« = -g 2 K{T + A VA )[^®iK 
I b 2 = - l -g 2 K{T + A VA ) 
F 2 = l -g 2 K(T + Asp) 



lu®lu 



ii®i„ 






C F 101 



(60) 
(61) 
(62) 



The total contribution including those from Fig. la', lb' and Id becomes 



2(I a 2 +I b 2 + r 2 ) = 2-g 2 K 



lu®lv 



ct/3;7<5 



101 



ij;kl 



xC F [T + A SP ] 



(63) 



which leads to 



A, 



2 N 2 -1 , m 
1 + g 2 — —{<T> + <A SP >} 



Ai 0) . 



(64) 



We again find that the vertex correction is multiplicative up to the one-loop level as opposed 
to the Wilson case 0. 

The contribution from the fermion self-energy has already been evaluated ||[7| and the 
total lattice renormalization factor is now obtained: 
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Z^ = (l-w 2 ) 2 Z 2 w Z 2 2 V r , 



(65) 



where 



T 



^r = l + 



16tt 2 

_£_ 

16n 2 



C h 



log(Aa) 2 + Ei 



— #rlog(Aa) 2 + vr 



v+ = 16n ^ tl [<T>(N + 1)+ « A VA » (N + 2)- « A SP »] 



V- = 

Vi = 

v 2 = 
with 



N 
16tt 2 (A^ + 1; 

N 
16tt 2 

16n 2 (N 2 - 1 



[< T > (N - 1)+ « A v A » (N - 2)+ « A S p >>] 



< T > (iV 2 - 1)+ << Ay A » N 2 - « A SP » + 5i logTr 2 , 
[<T> + « A SP »} +5 2 log7r 2 



N 



{ (N -l)(N -2) 



N 




(N + 1)(N + 


2) 


N 




(N + 2)(N- 


2) 


N 




A(N + 1)(N ■ 


-1) 



+ 



N 
The infrared singularity of < A x > is subtracted as 

d 4 k 



« A x » 



(2tt)< 



^(A;)^^) - c X7 - -;9(7r 2 - k 



(A' 



(66) 
(67) 

5 + logTr 2 , (68) 

5_logvr 2 , (69) 

(70) 

(71) 



with csp = 4 and cyA = 1- 

Numerical values of vr are evaluated by two independent methods. In one method the 
momentum integration is performed by a mode sum for a periodic box of a size L 4 after 
transforming the momentum variable through k^ = q^ — sing^. We employ the size L = 64 
for integrals. In the other method the momentum integration is carried out by the the Monte 
Carlo integration routine VEGAS, using 20 samples of 1000000 points each. We find that 
both results agree very well. Numerical values of v r are presented in Table [TT] as a function 
of M. 

We have to also calculate the corresponding continuum wave-function renormalization 
factor and vertex corrections in the MS scheme employing the same gauge and the same 
infrared regulator as the lattice case. For the present calculation it seems preferable to choose 



Dimensional Reduction(DRED) as the ultraviolet regularization, in which the loop momenta 
of the Feynman integrals are defined in D < 4 dimensions while keeping the Dirac matrices 
in four dimensions. In the DRED scheme we can use the same calculational techniques for 
the vertex corrections as the lattice case thanks to applicability of the Fierz transformation 
for the Dirac matrices. For the wave-function renormalization factor a simple calculation 
gives 



7MS 



1 + 



<T 



16vr 2 



c f 



log(A//i) 2 - 1/2 



where \x is a renormalization scale, 
corrections we obtain 



This result leads to T,{ 



MS 



(72) 



— 1/2. For the vertex 



V, 



MS 



16tt 2 



log(A//x) 2 + 1 



(73) 



giving v^ s 



dr- Here we should remark that the one-loop vertex corrections yield the 
evanescent operators which vanish in D = 4 for the DRED scheme ||. It is meaningless to 
give results without mentioning the definition of evanescent operators, because the constant 
terms at the one-loop level depend on the definition of the evanescent operators. Our choice 
is as follows: 



E. 



DRED 

± 



E\ 



DRED 



<Vfr(! - 7s) ® 7,(1 - 7s) - f 7m( 1 - 7s) ® 7/X 1 ~ 7s) 
Vy^ 1 - 7s) ® 7,(1 + 75) - f 7m(1 - 75) <8> 7^(1 + 75) 



(74) 
(75) 



where 5^ v is the D-dimensional metric tensor which emerges inevitably in the evaluation of 
the Feynman integrals. 

Combining these results with the previous lattice ones we obtain 



MS 



0*) 



'i-<) 2 zi 



Z r {na)O l f{l/a 



(76) 



where 



Z T (fj,a) 



[zrf 



\/ r Ms 



[z 2 yv r 



1 + 



16tt 2 



(&r - 2C F ) log(fia) 2 + z r 



Zt 



,MS 



Wr + 2C F {Sf b -E 1 }. 



(77) 
(78) 



Numerical values of zr are given in Table [TlT] and the results for the mean-field improved 



one, Zy , are also given in Table [TV . 



Although the results for the DRED scheme are presented here, it is an easy task to 
obtain those for the Naive Dimensional Regularization(NDR) scheme. In Appendix B we 
summarize the finite parts of the wave-function renormalization factor and vertex corrections 
in the NDR scheme. 
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IV. RENORMALIZATION FACTORS FOR THREE-QUARK OPERATORS 

The three-quark operators relevant to the proton decay amplitude are given by 

(O pd ) s = e abc ((q?) a T x ( q2 ) b ) (T Y (q 3 y) 5 , (79) 

where q c = —q T C^ 1 with C = 7072 is a charge conjugated field of q and Tx <S> Ty = 
Pr ® Pr, Pr ® Pl, Pl <8> Pri Pl <£> Pl- The summation over repeated color indices a, b, c is 
assumed. We should note that the domain-wall fermion action ([IJ) is transformed identically 
into that with the conjugated field by using transpose and matrix C. The resultant action 
and Feynman rules for the conjugated field is obtained by the replacement that 

igT A -> -ig(T A ) T , (80) 

where the superscript T means the transposed matrix. 

In order to evaluate the vertex corrections we consider the following Green function: 

((o PD ) 5 )f 01 = ((o PD ) s (rf )' a («,)j(«,);>, (si) 

where a, /3,7 and i,j,k are spinor and color indices respectively. Truncating the external 
quark propagators of (0 P d) we obtain the vertex function 

^■fc = n _„^3/ A (o) , A (i)V ife 



(1 " <f (A PD )T^ = (1 - ^oY (W + A^j Q ; ;57 , (82) 

where the trivial factor (1 — w^) 3 is factored out for the convenience. We suppress the 
external momenta pi since the renormalization factor does not depend on them. 
At the tree level the vertex function takes the form 

A?) = e^ k [V x ® IY] Q/3;57 , (83) 

where [T x <S> r y ] a/3;(57 = (r^^IY)^. 

The one-loop vertex corrections are shown in Figs. 2a, 2b and 2c, the sum of which gives 
the one-loop level vertex function 

A?t = £ -0ji (l a PD + I PD + I PD ) ■ (84) 

Using the notations in eqs.(^Sj), (|3"0D and fl3"2] ) the integrands Ipjf are written as follows: 



ra abkl rpA\T rpB 

1 PD — b \ 1 )ia 1 bj 



x {v,(k)G(k) + V,(k)G(k)}v x {G{k)V u {k) + G{k)V v {k)} <g> Y Y G^{k), (85) 



rb ibcrpArpB 

1 PD — e 1 bj 1 ck 



xT x {G(k)V u (k) + G{k)V v {k)} ® T Y [G(-k)V v {-k) + G(-A;)K(-A;)} G A *{k), (86) 

rC c a 3 c ( rpA\T rpB 

1 PD — e \ 1 )ia 1 ck 

x {V„(k)G(k) + V /I (k)G(k)}v x ® T y {£(fc)F„(fc) + G(fc)K(fc)} G^f (fc). (87) 



11 



A little algebra yields 

,N+1 



Ipd = g 2 

Ipd = 9 

I C PD = £ 



2N 
2 N + 1 

2N 
• iV + 1 ,,» 



e ijk K[T + A SP ][T x ^T 



Y\ 



: ijk K 



T(T X ® r y ) + cos 2 (k fl /2) sin 2 k^Txlal^) ® {Trials 



2N 



e ijk K T{T X <g> r y ) + cos 2 (A; M /2) sin 2 ^(t^x) ® (r y7o7/ ,; 



(88) 
(89) 
(90) 



where K, T and A S p are given in eqs. (|38D , (|39"D and ([51]). It is noted that a sum of I PD and 



/on becomes 



l PD 



2 N + 1 ijkT ^ 
g — e J if 



4 



2A 

2^+1 

2 i 

2A 



•yfe^" 



2T(T X <%> T Y ) + cos 2 (k^/2) sin 2 /c a ( 7Al7a r x + r x7a7A1 ) <8> (r y7ce7/U ) 
2T(r x (g) r y ) + cos 2 (£; M /2) sin 2 k a ({^ la + j a ^}T x ) ® (r y7a7/ J (91) 



for Tx = -Pr or Pj,, therefore no Fierz transformation is necessary to simplify the spinor 
structure of the total contribution. Finally we obtain 



Ipd + Ipd + 1 pd = 9 2 ^T £tlkK l T x ® T y ] [3T + A SP + 2A VA \ 



(92) 



Compared with the tree level result of eq.(83) we find that the vertex correction is multi- 
plicative up to the one-loop level: 



A 



PD 



1 + 9 2 ^~ {3 < T > + < A SP > +2 < A VA >} 



A 



(0) 
PD-- 



(93) 



where < X > (X = T, Ay a, Asp) are defined in eq.(|54|). We remark that in the Wilson case 
Opd mixes with other operators which have different chiral structures under renormalization 

Taking account of the contribution of the wave function the lattice renormalization factor 
for Opd is expressed as 



z% = (i - wlf'zTzTv, 



PD, 



where 



V, 



PD 



1 + 



16tt 2 



-5 PD log(Aa) 2 + v PD 



(94) 



(95) 



r rv = ^ i, r [3 < T > + « A SP » +2 « A VA »] + 5 PD logvr 2 (96) 



with Spo 



2N 
6(A + 1) 



. Numerical values for vpd, evaluated as before, are given in Table |TJ 



2A 
as a function of M. 

The corresponding continuum renormalization factors in the MS scheme are calculated 

employing the DRED scheme as the regularization in the Feynman gauge with the fictitious 

gluon mass A. The vertex correction for Opd is 
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V MS - 1 + 

V PD — L * 



16vr 2 



:S 



PD 



log(A//z) 2 + 1 



(97) 



giving t>p^ = Spd- We remark that in this case the evanescent operator does not appear at 
the one-loop level. 

Combining this result with the previous lattice one we finally obtain the relation between 
the operators ti% and 0%: 



OTM 



i 



;i - wifi^z, 



3/2 



Z PD (iia)Of D {l/a) 



(98) 



with 



(Z 2 MS)3/VMS 

ZMfia) ~ (Z 2 )^V PD 



1 + 



16vr 2 



(S PD - 3C F /2) log(/ia) 2 + z PD 



zp D = vfl-v PD -r°-C F {^ -EJ. 



(99) 
(100) 



We present numerical values for zp^ in Table [T[| and those for the mean-field improved one, 
z%, in Table |V. 



V. RENORMALIZATION FACTOR FOR B K 

As an application of results in the previous sections, we estimate a renormalization factor 
for the kaon B parameter Bk, defined by 



B 



K 



(TC\O + \K ) 



l(K v \A 4 \0)(0\A 4 \K°) 



101) 



with q 1 = q 3 = s and q 2 = q$ = d in + . 

Denoting the renormalization factor between the continuum Bk at scale // and the lattice 
one at scale l/o as Z BK (fia), we obtain 



Z BK (pia) 



[l-w^Z^Z+^a) __ Z+(fia) 



1 - u> ) 2 Z~ 2 Z A (fj,a) 2 Z A (/jay 



where 



Z+(jjui) = 1 + ^— [-41og(/xa) + 24 

l07T z 



from eq. (|77| ) in this paper, and 



Z A (fj,a) = 1 + 



16tt 2 



-z.a 



(102) 



(103) 



(104) 



from Ref. 0, so that 
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Z Bk {jm) = 1 + -f-^ [-41og(/ia) + ^ + - 2C F z A ] . (105) 

Note that 2U in Ref. @ is evaluated in the NDR scheme while the DRED scheme is used 
for z + in this paper. From the result in Appendix B we have 

za(DRED) = za(NDR) + 1/2, 2+(DRED) = ^+(NDR) + 3. (106) 

In Ref. || B K has been evaluated at (3 = 5.85, 6.0 with M = 1.7 and (3 = 6.3 with 
M = 1.5, using domain-wall QCD with the quenched approximation. Here we explicitly 
calculate Zg K (i/,a) for these parameters. From Table [Til] and the previous result 0, z + = 
-41.854(-42.399), z A = -17.039(-16.827) and z + -2C F z A = 3.583(2.473) for M = 1.7(1.5) 
in the DRED scheme, and z + = -44.854(-45.399), z A = -17.539(-17.327) and z + - 
2C F z A = 1.917(0.8063) for M = 1.7(1.5) in the NDR scheme. Taking \x = l/o and g 2 = 
gjj^(l/a), estimated by the formula 

-5— (l/o) =P-- 0.13486 (107) 

9jls 6 

for the quenched QCD with P being the average value of the plaquette, we have Zb k = 
1.053 (1.029), 1.049 (1.026) and 1.030 (1.010) at (3 = 5.85, 6.0 and 6.3, respectively, in the 
DRED (NDR) scheme. Sizes of one-loop corrections for B^ are not so large, 1 — 5%, at these 
j3 values even without mean-field improvement, since the large contribution, which comes 
from a (1 — wo)Z w factor, cancels out in the ratio of + and A\. 

If we employ the mean-field improvement by replacing M —>■ M = M + 4(u — I) with 
u = P 1 /*, we obtain Z Bk = 1.018 (0.994), 1.017 (0.994) and 1.009 (0.988) at (3 = 5.85, 6.0 
and 6.3, respectively, in the DRED (NDR) scheme. See Appendix A for some remarks. 

Note that there is no mean-field improvement factor for Bk in actual simulations since, 
as mentioned before, it is defined by the ratio. Therefore the difference between values 
of Zb k with and without mean-field improvement comes from higher order ambiguity in 
perturbation theory. 

Necessary informations for the analysis in this section are given in Table |V|, together 
with values of Z Bk . 

VI. CONCLUSION 

In this paper we have calculated the one-loop contributions for the renormalization fac- 
tors of the three- and four-quark operators in DWQCD. We have demonstrated that the 
three- and four-quark operators in DWQCD can be renormalized without any operator mix- 
ing between different chiralities as opposed to the Wilson case. This desirable property in 
DWQCD would practically surpass the cost of the introduction of an unphysical fifth dimen- 
sion. The numerical values for the finite parts zx with X = ±, 1, 2, PD settle in reasonable 
magnitude with the mean-field improvement, while unimproved values are rather large in 
general. 

In this work we do not treat the operators which yield the so-called "penguin" diagram. 
It seems feasible to carry out the calculation of their renormalization factors, which we leave 
to future investigation. 
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APPENDIX A: MEAN-FIELD IMPROVEMENT 



The mean-field improvement [11 j in our paper uses 



u 



1 - 



9 2 C, 



T 



with T = 0.15493, which is the value for the link in Feynman gauge. It may be better to 
use u from K c or plaquette in DWQCD. In that case 



u 



1 



9 2 C, 



(T + ST), 



where ST = 0.00793 for K c , or ST = —0.02993 for plaquette. Accordingly we have to modify 
renormalization factors as follows: 



2w 



Z M> (T + ST) = Z M* (T) + 

z% IF (T + ST) = z% IF (T) + 16tt 2 x ST/2 



16tt 2 x 25T, 



^r, bilinear 



(T + ST) = z- 



r,iv me ,r(T) + l6^x5T/2, 

MF 



4\Uu,AT + ST) = z^_ quark (T) + 16^ x ST/2 x 2C F , 
4T^AT + ST) = z^l quark (T) + 167T 2 x ST/2 x *-C F . 



(108) 

(109) 
(110) 
(111) 

(112) 



APPENDIX B: NAIVE DIMENSIONAL REGULARIZATION(NDR) 

In this appendix we compile the finite part of the renormalization constant in the MS 
subtraction scheme with the Naive Dimensional Regularization: 



Ef B = 1/2, 



,MS 



.MS 



MS 



, + x{3/2 _^ } , 



cL x {3/2 
' 5 1 x 1 3/2 



N -2 
2iV-3 



}• 



N + 2 
2(N 2 - 5) 



N 2 



\ 



N 



x 3, 



v$B = S PD x 2/3, 



S 2 x 3/4 
S 2 x 1/2 



(113) 
(114) 

(115) 

(116) 
(117) 
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where Vij with i,j = 1, 2 is a matrix, which represents the mixing of the finite part for 0\^- 
The one-loop vertex corrections for Or (T = ±, 1, 2) require to specify their evanescent 
operators, which originates from the property that the Fierz transformation can not be 
defined in the NDR scheme. We employ 

E™ = 7 p 757m(1 " 7s) ® 7m(1 ~ 7 5 )7<57p - (2 - D) 2 lfl (l - 7s ) ® 7 „(1 - 7s ), (118) 

^i N 2 m = 7p757m(1 - 7s) ® 7m(1 + ls)ltl P - D 2 ltl (l - 7s ) <g> 7m (1 + 7s ), (119) 

where -D is the reduced space-time dimension. On the other hand, the evanescent operator 
does not appear in the one-loop vertex correction of Opr>. 

For later convenience, values of the finite part of quark bilinear operators are also given 
here. For NDR scheme 

^1 = 0, 4? = 5/2, ^=1/2, (120) 

while for DRED scheme 

*S = l/2, 4? = 7/2, ^="1/2, (121) 

where the evanescent operators are 

< RED = <W - j7„ (122) 



D 



<75 = Vvfi - T 7,75- (123) 
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TABLES 



TABLE I. Color factors for I^' b ' c (L = ±, 1, 2). 



-n T AB jAB jAB 

^ ' J _a J b ,J _c 

± T A T B ®\ ± T A QT B T a ®T b ± T A QT B T a ®T b ± T^T^Ql 

1 -NT A T B ®1 + T A QT B -NT a ®T b + T a QT b -NT a ®T b + T A T B Ql 

2 T A QT B T A QT B T A T B Ql 



TABLE II. Numerical values for V r (T = ±, 1, 2, PD) as a function of M. 



M 


v+ 


V_ 


V! 


v 2 


V PD 


0.05 


13.9096(8) 


10.847(8) 


13.3992(19) 


8.805(12) 


8.646(5) 


0.10 


13.5696 


11.537 


13.2309 


10.182 


8.992 


0.15 


13.2941 


12.098 


13.0948 


11.301 


9.273 


0.20 


13.0548 


12.587 


12.9768 


12.275 


9.518 


0.25 


12.8391 


13.029 


12.8708 


13.155 


9.740 


0.30 


12.6404 


13.438 


12.7734 


13.970 


9.946 


0.35 


12.4542 


13.822 


12.6822 


14.734 


10.139 


0.40 


12.2775 


14.188 


12.5960 


15.462 


10.323 


0.45 


12.1083 


14.539 


12.5135 


16.160 


10.499 


0.50 


11.9449 


14.880 


12.4341 


16.837 


10.671 


0.55 


11.7861 


15.212 


12.3571 


17.496 


10.838 


0.60 


11.6307 


15.538 


12.2819 


18.143 


11.002 


0.65 


11.4779 


15.859 


12.2081 


18.780 


11.164 


0.70 


11.3269 


16.178 


12.1354 


19.412 


11.325 


0.75 


11.1770 


16.495 


12.0634 


20.041 


11.485 


0.80 


11.0275 


16.813 


11.9917 


20.670 


11.645 


0.85 


10.8779 


17.131 


11.9201 


21.300 


11.806 


0.90 


10.7276 


17.452 


11.8484 


21.935 


11.968 


0.95 


10.5760 


17.777 


11.7762 


22.578 


12.133 


1.00 


10.4225 


18.107 


11.7033 


23.230 


12.300 


1.05 


10.2659 


18.437 


11.6278 


23.885 


12.466 


1.10 


10.1076 


18.790 


11.5547 


24.579 


12.646 


1.15 


9.9443 


19.139 


11.4768 


25.269 


12.822 


1.20 


9.7768 


19.505 


11.3981 


25.990 


13.007 


1.25 


9.6037 


19.880 


11.3165 


26.732 


13.198 


1.30 


9.4244 


20.272 


11.2323 


27.504 


13.396 


1.35 


9.2375 


20.680 


11.1446 


28.309 


13.603 


1.40 


9.0419 


21.109 


11.0530 


29.153 


13.820 


1.45 


8.8361 


21.560 


10.9567 


30.042 


14.049 


1.50 


8.6183 


22.037 


10.8547 


30.983 


14.291 


1.55 


8.3863 


22.547 


10.7464 


31.987 


14.550 


1.60 


8.1375 


23.093 


10.6301 


33.063 


14.827 


1.65 


7.8685 


23.683 


10.5043 


34.227 


15.127 


1.70 


7.5747 


24.328 


10.3668 


35.496 


15.454 


1.75 


7.2502 


25.038 


10.2149 


36.897 


15.814 


1.80 


6.8864 


25.832 


10.0440 


38.463 


16.216 


1.85 


6.4706 


26.737 


9.8482 


40.247 


16.675 


1.90 


5.9812 


27.794 


9.6168 


42.337 


17.210 


1.95 


5.3749 


29.093 


9.3281 


44.907 


17.867 
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TABLE III. Numerical values for z-p (r = ±, 1, 2, PD) as a function of M. 



M 


z + 


z_ 


21 


Z2 


zpd 


0.05 


-49.908(10) 


-40.845(11) 


-48.397(8) 


-34.803(15) 


-32.144(7) 


0.10 


-49.332 


-41.300 


-47.994 


-35.945 


-32.314 


0.15 


-48.847 


-41.651 


-47.647 


-36.853 


-32.437 


0.20 


-48.416 


-41.949 


-47.338 


-37.637 


-32.539 


0.25 


-48.025 


-42.215 


-47.057 


-38.341 


-32.629 


0.30 


-47.663 


-42.461 


-46.796 


-38.993 


-32.713 


0.35 


-47.325 


-42.693 


-46.553 


-39.605 


-32.792 


0.40 


-47.007 


-42.918 


-46.326 


-40.191 


-32.870 


0.45 


-46.706 


-43.137 


-46.111 


-40.758 


-32.948 


0.50 


-46.419 


-43.354 


-45.908 


-41.311 


-33.027 


0.55 


-46.145 


-43.571 


-45.716 


-41.855 


-33.108 


0.60 


-45.883 


-43.790 


-45.534 


-42.395 


-33.191 


0.65 


-45.631 


-44.012 


-45.361 


-42.933 


-33.279 


0.70 


-45.388 


-44.239 


-45.196 


-43.473 


-33.370 


0.75 


-45.153 


-44.472 


-45.040 


-44.017 


-33.467 


0.80 


-44.927 


-44.712 


-44.891 


-44.569 


-33.570 


0.85 


-44.708 


-44.961 


-44.750 


-45.130 


-33.679 


0.90 


-44.496 


-45.220 


-44.617 


-45.704 


-33.795 


0.95 


-44.290 


-45.491 


-44.490 


-46.292 


-33.918 


1.00 


-44.091 


-45.776 


-44.372 


-46.899 


-34.051 


1.05 


-43.898 


-46.070 


-44.260 


-47.517 


-34.190 


1.10 


-43.710 


-46.392 


-44.157 


-48.181 


-34.347 


1.15 


-43.528 


-46.723 


-44.061 


-48.853 


-34.510 


1.20 


-43.352 


-47.079 


-43.973 


-49.565 


-34.688 


1.25 


-43.180 


-47.457 


-43.893 


-50.308 


-34.880 


1.30 


-43.014 


-47.862 


-43.822 


-51.094 


-35.088 


1.35 


-42.853 


-48.296 


-43.760 


-51.924 


-35.315 


1.40 


-42.697 


-48.763 


-43.708 


-52.808 


-35.561 


1.45 


-42.545 


-49.269 


-43.666 


-53.751 


-35.831 


1.50 


-42.399 


-49.818 


-43.635 


-54.763 


-36.127 


1.55 


-42.257 


-50.417 


-43.617 


-55.857 


-36.453 


1.60 


-42.119 


-51.074 


-43.611 


-57.044 


-36.813 


1.65 


-41.985 


-51.800 


-43.621 


-58.343 


-37.214 


1.70 


-41.854 


-52.607 


-43.646 


-59.776 


-37.663 


1.75 


-41.725 


-53.513 


-43.690 


-61.372 


-38.170 


1.80 


-41.595 


-54.541 


-43.753 


-63.172 


-38.748 


1.85 


-41.460 


-55.727 


-43.838 


-65.237 


-39.417 


1.90 


-41.311 


-57.124 


-43.946 


-67.666 


-40.207 


1.95 


-41.121 


-58.840 


-44.074 


-70.652 


-41.177 
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TABLE IV. Numerical value for z^ F {V = ±, 1, 2, PD) as a function of M. 



M 


z + 


Z MF 


z l 


z 2 


~MF 
Z PD 


0.05 


-17.287(10) 


-8.224(11) 


-15.776(8) 


-2.182(15) 


-7.679(7) 


0.10 


-16.712 


-8.679 


-15.373 


-3.324 


-7.848 


0.15 


-16.226 


-9.030 


-15.027 


-4.232 


-7.972 


0.20 


-15.796 


-9.328 


-14.718 


-5.016 


-8.074 


0.25 


-15.404 


-9.594 


-14.436 


-5.720 


-8.164 


0.30 


-15.042 


-9.840 


-14.175 


-6.372 


-8.247 


0.35 


-14.705 


-10.073 


-13.933 


-6.98 


-8.326 


0.40 


-14.386 


-10.297 


-13.705 


-7.57 


-8.404 


0.45 


-14.085 


-10.516 


-13.490 


-8.13 


-8.482 


0.50 


-13.799 


-10.734 


-13.288 


-8.69 


-8.561 


0.55 


-13.525 


-10.951 


-13.096 


-9.23 


-8.642 


0.60 


-13.262 


-11.169 


-12.913 


-9.77 


-8.726 


0.65 


-13.010 


-11.391 


-12.740 


-10.3 


-8.813 


0.70 


-12.767 


-11.618 


-12.575 


-10.8 


-8.905 


0.75 


-12.532 


-11.851 


-12.419 


-11.3 


-9.002 


0.80 


-12.306 


-12.091 


-12.270 


-11.9 


-9.104 


0.85 


-12.087 


-12.340 


-12.129 


-12.5 


-9.213 


0.90 


-11.875 


-12.600 


-11.996 


-13.0 


-9.329 


0.95 


-11.670 


-12.871 


-11.870 


-13.6 


-9.453 


1.00 


-11.470 


-13.155 


-11.751 


-14.2 


-9.585 


1.05 


-11.278 


-13.449 


-11.639 


-14.8 


-9.725 


1.10 


-11.089 


-13.772 


-11.536 


-15.5 


-9.882 


1.15 


-10.908 


-14.103 


-11.440 


-16.2 


-10.044 


1.20 


-10.731 


-14.459 


-11.352 


-16.9 


-10.223 


1.25 


-10.559 


-14.836 


-11.272 


-17.6 


-10.414 


1.30 


-10.393 


-15.241 


-11.201 


-18.4 


-10.623 


1.35 


-10.232 


-15.675 


-11.139 


-19.3 


-10.849 


1.40 


-10.076 


-16.143 


-11.087 


-20.1 


-11.096 


1.45 


-9.925 


-16.648 


-11.045 


-21.13 


-11.365 


1.50 


-9.778 


-17.197 


-11.014 


-22.14 


-11.661 


1.55 


-9.636 


-17.796 


-10.996 


-23.23 


-11.987 


1.60 


-9.498 


-18.453 


-10.991 


-24.42 


-12.347 


1.65 


-9.364 


-19.179 


-11.000 


-25.72 


-12.749 


1.70 


-9.233 


-19.986 


-11.025 


-27.15 


-13.198 


1.75 


-9.104 


-20.892 


-11.069 


-28.75 


-13.705 


1.80 


-8.975 


-21.920 


-11.132 


-30.55 


-14.283 


1.85 


-8.840 


-23.106 


-11.217 


-32.61 


-14.952 


1.90 


-8.690 


-24.503 


-11.326 


-35.04 


-15.742 


1.95 


-8.500 


-26.219 


-11.453 


-38.03 


-16.711 
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TABLE V. Renormalization factor for Bk (1/a) at some parameters. 



p 


5.85 


6.0 


6.3 


M 


1.70 


1.70 


1.50 


P 


0.57506 


0.59374 


0.62246 


^(V«) 


2.3484 


2.1792 


1.9278 


u 


0.87082 


0.87781 


0.88823 


M 


1.20 


1.20 


1.05 




DRED NDR 


DRED NDR 


DRED NDR 


z+ 


-41.854 -44.854 


-41.854 -44.854 


-42.399 -45.399 


ZA 


-17.039 -17.539 


-17.039 -17.539 


-16.827 -17.327 


z + - 2Cfza 


3.583 1.917 


3.583 1.917 


2.473 0.806 


Z Bli {^a = 1) 


1.053 1.029 


1.049 1.026 


1.030 1.010 




-17.033 -20.033 


-17.033 -20.033 


-17.580 -20.580 


~MF 
Z A 


-6.853 -7.353 


-6.853 -7.353 


-6.864 -7.364 


zf F - 2C F zf F 


1.242 -0.425 


1.242 -0.425 


0.724 -0.943 


Zg K F (»a = 1) 


1.018 0.994 


1.017 0.994 


1.009 0.988 
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FIGURES 
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FIG. 1. One-loop vertex corrections for the four-quark operator. a,j3, 7,(5 and i,j,k,l label 
Dirac and color indices respectively. 
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FIG. 2. One-loop vertex corrections for the three-quark operator, a, /3, 7 and i,j,k label Dirac 
and color indices respectively. 
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